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Abstract 

We study a one-dimensional model for granular gases, the so-called Inelastic Maxwell 
Model. We show theoretically the existence of stationary solutions of the unforced 
case, that are characterized by an infinite average energy per particle. Moreover, 
we verify the quasi-stationarity of these states by performing numerical simulations 
with a finite number of particles, thereby highlighting truncated Levy distributions 
for the velocities. 
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1 Introduction 



High energy tails are generic for the velocity distributions of fluidized granular 
media, i.e. low density systems of macroscopic particles performing inelastic 
interactions. Contrary to the specific shape of the tail which depends on the 
details of the grains, overpopulation seems to be a general feature of inelastic 
systems, which expresses their deep non-equilibrium nature. Recently, classes 
of simplified models, the so-called Inelastic Maxwell Models (IMM) have been 
introduced. These models derive from a mathematical simplification of the 
non-linear Boltzmann equation for inelastic hard spheres, assuming that the 
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collision rate may be replaced by an average rate, independent of the relative 
velocity of the colliding particles. This model introduces mathematical sim- 
plifications which may help to highlight the origin of the anomalous velocity 
distributions and serve as a simple tool to study far-from-equilibrium systems. 
Quoting M.H. Ernst, one can consider that "what harmonic oscillators are for 
quantum mechanics and dumb-bells for polymer physics, is what elastic and 
inelastic Maxwell models are for kinetic theory," (ref.(0)). 

Usually, people focus on two kinds of solutions of the inelastic Boltzmann 
equation, namely scaling solutions describing the approach towards total rest 
state and heated stationary solutions. In the following, however, we focus on 
another class of solutions of inelastic gases, which are stationary states of the 
unforced case. Before going further, let us point the recent exhaustive study 
of Ben-Nairn and Machta (ref.(2)), that focuses independently on the same 
issue. Our analysis is done in the context of the one-dimensional IMM, where 
we show that stable Levy distributions play a central role. The reasons for 
studying these stationary Levy velocity distributions are multiple. First, they 
correspond to exact solutions for an inelastic gas, and this fact is very rare. By 
discussing the physical relevance of these states, whose mean energy diverges, 
we show that there exists physical solutions which are very close to these 
stationary states and which may be considered quasi stationary. This paper is 
also motivated by the revived interest regarding generalized thermodynamical 
equilibrium theories. Indeed, long tail distributions have been found in diverse 
fields of science and are the object of a growing number of researches. From the 
point of view of statistical physics, several attempts have been made in order 
to show physical mechanisms underlying these anomalous behaviours, such 
as Tsallis's non-extensive entropy-based statistics (ref. (0)). The study of the 
above Levy distributions could help to build a thermodynamic-like foundation 
to processes such as anomalous diffusion, fractal dynamics... This problem, 
which has first been questioned by Montroll and Schlesinger (ref. (jj)), has 
recently lead to study simple kinetic models in order to show equilibrium-like 
features of Levy distributions (ref. (0), (0)). In this work, we will also discuss 
the range of validity of such Levy states, and show that they are limited to 
the context of Maxwell-like kinetic models. 



2 Levy velocity distribution 



The one- dimensional IMM has been first introduced by Ben-Nairn and Krapivsky 
(ref.(0)), and derives from a simplification of the Boltzmann equation for in- 
elastic hard rods. The collision rule between grains 1 and 2 reads: 

(1 + a) , (1 + a) /iN 

v 1 = v 1 v 12 v 2 = v 2 H v 12 (1) 

la la 
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where the prime velocities are the pre-collisional velocities. The inelasticity is 
accounted through the so-called normal restitution coefficient a restricted to 
the interval [0:1]. Standard procedures, such as molecular chaos hypothesis, 
lead to the one- dimensional inelastic Boltzmann equation. The IMM follows 
after the replacement of the collision rate by an average rate and a rescaling 
of the time scale: 

()Iiri: t) +f(v 1 ;t)= J dv 2 (-)f(v[;t)f(v' 2 ;t) (2) 



dt J a 



where one has furthermore assumed that the system is and remain homo- 
geneous in the course of time. Finite energy solutions of (2) have a natural 
tendency to reach the asymptotic solution S(v), expressing the natural cool- 
ing down due to the dissipative interactions. These solutions have been proven 
to exhibit multiscaling (ref.(0)), i.e the system is characterized by an infinite 
set of independent cooling rates. Moreover, as shown by Baldassarri et al. 
(ref.(ji|)), multiscaling leads to a power law scaling velocity solution. 

It is well known that the Maxwell like collision operator of (2) is very similar to 
a convolution process (ref.g)). This property is clear after using the Bobylev 
Fourier transform method to the kinetic equation: 

^ + ^( M ) = «(i±^) MM (^i) M ) (3) 



where 4>(k; t) is the characteristic function of the velocity distribution, defined 
by: </)(k; t) — j- J dvf(v,t)e lkv . This expression reveals the mechanisms which 
dominate the long time dynamics of the kinetic model. Indeed, the evolution of 
<f)(k*; t) is determined by the values of <p(k; t) for \k\ < \k*\ (this comes directly 
from the inequalities < \k*\ and \(^r)k*\ < \k*\). Consequently, as 

time increases, <p(k*; t) is determined by the initial condition (p(k; 0) for smaller 
and smaller values of |fc| and, in the long time limit, the characteristic functions 
are dominated by their initial behaviour at k — > 0. If the initial energy diverges, 
and if the small k development writes 4>(k) = l + a(t)\k\, one easily shows that 
d t a(t) = 0, i.e. the quantity a describing the small k singularity of 4>{k) is a 
constant of motion. In that case, the stationary solution of (3) is the Levy 
distribution: 

(j)(k) = e- a|fe| (4) 



which, in velocity space rewrites: 

M = ~ ( 2 I 2^ (5) 

ii (a 2 + v 2 ) 
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Let us stress that this class of states is characterized by an infinite energy, a 
condition which is a priori necessary in order to make compatible the inelastic- 
ity of the interactions and the stationarity of the velocity distribution. More- 
over, analogies with the Central-Limit theorem suggest that (5) is the asymp- 
totic stationary solution for a large class of initial distributions, thereby draw- 
ing a parallel with the Maxwell-Boltzmann equilibrium distribution (ref . (fioh ) . 
One should also remark that the quantity a also shares some properties with 
an equilibrium temperature, such as the equipartition property in the case of 
mixtures (ref.fjH)). 



3 Physical relevance and simulations of Levy distributions 



In this section, we would like to discuss the physical relevance of the station- 
ary Levy velocity distributions derived in the previous section. At first sight, 
the only physical stationary solution should be the total rest state, <p(k) = 1 
(e -a ' fc !, with a=0), because an infinite average energy for the velocity distribu- 
tion sounds unphysical. It is nonetheless interesting to point out that, despite 
this infinite average behaviour, each realization of the statistical process is 
finite, i.e. each particle carries a finite but unbounded energy. Moreover, phys- 
ical systems are composed by a finite number of constituents. Combining these 
two remarks, one expects that the average energy of a physical system is finite 

l-a 2 

and decreases because of dissipation according to Haff law: T(t) = Tq e ~ *, 
even if its velocities are initially Levy distributed. In order to understand this 
paradox, we have performed computer simulations of the Inelastic Maxwell 
Model. This can be achieved by using the usual stochastic interpretation of 
the Boltzmann equation (ref.()13»). To do so, we initially distribute velocities 
according to a Lorentz distribution, and check that the energy of the system 
decreases like the exponential law prescribed by theory. The velocity distribu- 
tion function is also shown to converge towards the Dirac delta distribution in 
the long time limit, as expected. In contrast, a careful study of the evolution 
of the core of the velocity distribution shows that it remains stationary during 
a non-negligible time interval (fig. 1). This property, together with the finite- 
ness of the number of particles, allows one to make a link between the IMM 
and truncated Levy distributions (TLD). This family of distributions has been 
introduced by Mantegna and Stanley (ref. fllh. flZ)) in order to justify the 
use of Levy-like distributions in finite-size systems. They are defined by: 



ftrunc( x ) — fh{ x ) if \ x \ < I 

ftrunc( x )=Q if \ x \ > I (6) 

where 1 is the cutoff length and /f is the Levy distribution of order (3 defined to 
be the inverse Fourier transform of e~ a ^ P . In a simulation of the IMM, the cut- 
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Fig. 1. Evolution of the core of the velocity distribution in log scale. The cutoff 
length of the Levy distribution is 10000. The inelasticity is a = 0.5 and the number 
of particles 10000000. We plot the velocity distribution after 1, 8 and 12 collisions 
per particle. 

off length is naturally imposed by the upper boundary value of the energy of 
the particles. In order to evaluate the stationarity of the Levy core distribution, 
we have measured the evolution of /(0; t). Asymptotically, this variable should 

behave like f(0\t) ~ e^^*, which is predicted from the scaling solution of 
Baldassari (ref.(J8|)). In the simulations, we observe the transition between two 
linear regimes in logarithmic scale (fig. 2): an initial regime during which /(0; t) 
is constant which corresponds to the quasi stationary regime of the core of the 
velocity distribution, and a linearly increasing regime, which comes from the 
energy decrease in the system. Amazingly, while f(0;t) increased of only ten 
percent during ten mean collision times, the energy decreased by a factor 100 
during the same time period. 



4 Conclusions 



As a conclusion, we would like to discuss the range of validity of the sta- 
tionary velocity distributions described above. Indeed, the Inelastic Maxwell 
Model studied in this paper is exceptional because of its strong similarities 
with random walk theory, which do not apply in general to more complex 
kinetic models. Therefore, it is important to recall a qualitative difference be- 
tween elastic and inelastic systems, namely the fact that the asymptotic solu- 
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Fig. 2. Evolution of the probability of return f(0;t) for the same system as in 
figure (1). The lines correspond to the theoretical predictions of the two regimes. 
The asymptotic saturation is due to the finite size of the measurement cells. 



tions of inelastic gases depend on their modeling while elastic systems always 
converge towards the Maxwell-Boltzmann distribution. This universality relies 
on the conservation of energy and on the time symmetric dynamics of elastic 
systems. Altogether, these properties imply the emergence of detailed balance 
at equilibrium. In contrast, it is straightforward to prove that stationary Levy 
distributions do not exhibit detailed balance (ref.([lfi)). thereby showing the 
dependence of the stationary solution on the collisional kernel. This result 
confirms that Levy distributions are associated to a summation of random 
variables, while the Gaussian also arises from more general principles. 

Nevertheless, Maxwell models constitute an ideal framework, due to their 
mathematical simplicity, in order to study non-trivial features of the Boltz- 
mann equation, such as non-equipartition of ener gy i n inelastic mixtures, or 
infinite energy solutions of the elastic case (ref. ()14 )). Moreover, they lead 
to a natural derivation of fractional Fokker-Planck equations (ref. (0)), and 
to the emergence of Levy distributions, whose interest in physics is therefore 
enhanced. In this paper, we also discussed the physical relevance of such infi- 
nite mean energy solutions by performing computer simulations of the model. 
Given the finiteness of the number of particles in the system, we showed that 
these solutions are equivalent to truncated Levy velocity distributions, i.e. to 
a finite mean energy solution. It is important to stress that such a solution 
emerges from the dynamics itself, and not from an ad hoc assumption for 
the truncature. Moreover, computer simulations indicate that, while the total 
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energy decays exponentially fast, as predicted by theory, the core of the distri- 
bution remains stationary during an arbitrary long period of time. Finally, we 
would like to point the recent work by Ben-Nairn and Machta (ref.(Q)), who 
propose an experimental realization of infinite energy steady states by inject- 
ing energy at large velocities. This experiment should open new perspectives 
in kinetic theory, and motivate detailed studies of this family of solutions. 

We would like to thank Dr. Figueiredo for fruitful discussions. R.L. would also 
like to thank especially Dr. Ausloos and Dr. Mareschal for invitations at the 
Universite de Liege and at CECAM respectively. This work has been done 
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References 

[I] M.H. Ernst and R. Brito, Europhys. Lett, 58, 2002, 182 

[2] E. Ben-Nairn and J. Machta, Phys. Rev. Lett., 94, 2005, 138001 
[3] C. Tsallis, J. Stat. Phys., 52, 1988, 479 

[4] E.W. Montroll and M.F. Schlesinger, J. Stat. Phys., 32, 1983, 209 
[5] E. Barkai, Phys. Rev. E, 68, 2003, 055104(R) 

[6] D.H. Zanette and M.H. Montemurro, Phys. Lett. A, 316, 2003, 184. 
[7] E. Ben-Nairn and P.L. Krapivsky, Phys. Rev. E, 61, 2000, R5 
[8] A. Baldassarri, U.M.B. Marconi and A. Puglisi, Europhys. Lett, 58, 2002, 
14 

[9] M.H. Ernst, Phys. Rep., 78, 1981, 1. 

[10] R. Lambiotte, phD thesis, Universite Libre de Bruxelles. 

[II] R.N. Mantegna and H.E. Stanley, Phys. Rev. Letters, 73, 1994, 2946 
[12] M.F. Schlesinger, Phys. Rev. Lett, 74, 1995, 4959 

[13] G.A. Bird , Molecular Gas Dynamics, Clarendon Press, 1976 
[14] A.V. Bobylev and C. Cercignani, J. Stat. Phys., 110, 2003, 333 



7 



